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Revisiting numerical real-space renormalization group for quantum lattice systems 
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Although substantial progress has been achieved in solving quantum impurity problems, the 
numerical renormalization group (NRG) method generally performs poorly when applied to quantum 
lattice systems in a real-space blocking form. The approach was thought to be unpromising for most 
lattice systems owing to its flaw in dealing with the boundaries of the block. Here the discovery 
of intrinsic prescriptions to cure interblock interactions is reported which clears up the boundary 
obstacle and is expected to reopen the application of NRG to quantum lattice systems. While the 
resulting RG transformation turns out to be strict in the thermodynamic limit, benchmark tests of 
the algorithm on a one-dimensional Heisenberg antiferromagnet and a two-dimensional tight-binding 
model demonstrate its numerical efficiency in resolving low-energy spectra for the lattice systems. 
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Since the success in solving the Kondo problem by 
Wilson Q, there had been lots of attempts in apply¬ 
ing the numerical renormalization group (NRG) method 
to treat other quantum many-body problems in a simi¬ 
lar way. However, the numerical RG algorithm based on 
real-space schemes met its Waterloo as it performed very 
poorly in subsequent several applications. In particular, 
the approach was shown to give inaccurate results for 
most interacting lattice systems n and for the Ander¬ 
son localization problem SI- The difficulty, identified 
by White and Noack [?] via a simple one-dimensional 
tight-binding model, lies in that the approach is flawed 
in its treatment of boundaries of a block. 

The inability to apply the NRG method to quan¬ 
tum lattice systems is a heavy loss to the research of 
condensed matter physics. Although a closely related 
method, the density matrix RG algorithm j8HlO(. was 
proposed later and has been shown very effective in 
achieving the ground-state energy for interacting lattice 
systems, there exist strong restrictions of this method in 
its application to excitation spectra or to systems with 
high spatial dimensions. The aim of this letter is to show 
that the flaw of the NRG method exposed previously 
can be eliminated via intrinsic prescriptions to cure in¬ 
terblock interactions and the derived “regularized” real- 
space blocking version of the NRG scheme is able to yield 
reliable results for quantum lattice systems. 

In the standard real-space version of the NRG ap¬ 
proach, one starts from a block Hamiltonian Hl of L sites 
and diagonalizes it exactly. By keeping a certain amount 
of the lowest-energy states {| ftj ), j = 1, • • •, to}, one then 
uses them to construct a Hamiltonian for a larger system 
composed of two such blocks: H 2 l = + v\ 2 + h\ 2 } , 

where v\ 2 denotes the interblock coupling. The primi¬ 
tive algorithm through projecting H 2 l simply onto to 2 
tensor product states \ipj) <8> | ipy) fails to achieve reli¬ 
able results. The new version of the NRG scheme here 
adopts a slightly different route: one should start by di¬ 
agonalizing a pair of block Hamiltonians, e.g., the one 
Hl with an open boundary condition and the one with 


a periodic form, H V L . Let us denote correspondingly the 
low-lying eigenstates of H £ by (| V'p}- The effective per¬ 
formance of the NRG algorithm resides in that, instead 
of using the set of tensor product states, the new pair of 
compound Hamiltonians, the open H 2 l and the periodic 
H!lj = H 2 l + V 21 , are constructed and diagonalized in 
virtue of the following two sets of to 2 states 


1 ^ jj') ~ U a dj '■ \lpj) <S> 

(1) 


(2) 


respectively. The transformation U a dj here is the key 
prescription introduced to cure the interblock coupling. 
It is formally a range-L operator acting on adjacent seg¬ 
ments of the compound block, namely, on the sites from 
fL + 1 to |L to generate |jy), and on the regions of 
the above sites and the sites from ~L + 1 to ^L sepa¬ 
rately so as to generate |'1/G,). As will be elucidated in 

detail, U a dj involves two different expressions, and 
U^dj'i the latter takes the form of V^ IV’JXVjf l> in 

which \il>f) = denotes a state through translat¬ 

ing | ipj) by L/2 sites and T\il) = |*i) has been assumed. 

The theoretical foundation of the above formulated 
prescription rests on the following analysis about the ut¬ 
most case in which the size L of the block is sufficiently 
large so that the system Hl can be viewed as thermally 
extensive. By dividing the block into two subsystems 
with intermediate coupling, Hl = hi + Vi r + h ri one can 
make an assertion that the low-lying spectra of Hl and 
of the disconnected Hamiltonian hi + h r are identical: 
E n = E' ni + E ^ 2 , where E l ni and E^ 2 denote the spectra 
of the left half hi and of the right half h r , respectively. 
This is simply understood since the partition function 
of the block, Z = Tre~ l3HL , should fulfill Z = Z/Z,. for 
any finite temperature 1//I; the latter identity, in which 
Z= Tre~ /3hl ’ r , actually accounts for the additivity of 
thermodynamic potentials of macroscopic thermal sys¬ 
tems. As Hl and hi + h r possess the same low-lying 


2 


spectra, the low-energy portion of the two systems could 
be linked by an isometry map, which turns out to de¬ 
lineate the first form of the prescribed transformation: 
ujutj = EJV’nX^n|, With I ip n ) = \ip ni n 2 ) denoting the 
low-lying eigenstates of hi + h r . That is, the low-energy 
portion of Hl could be obtained via U^(hi + h r )\U^]\ 
and the effect of the intermediate coupling vi r hence is 
equally described by the transformation U^. 

It is then crucial to note that the formally range-L 
indeed has an effective range less than L as the block size 
L increases asymptotically. Especially, for the infinite Hl 
with a gapped ground state, the range of turns out to 
be local since the evolution operator generated by Hl in 
a finite period r, formed as e~ lhlT ®e~ lhrT sandwiched by 
U^2 and could be simulated by a quantum circuit 

11] with finite depth. The latter representative, which 
describes the time evolution operator in virtue of a few 
layers of piecewise local unitary operators, has also been 
exploited to characterize quantum phases and topological 
orders for quantum many-body systems fl2l . [l3| . 


Consequently, one can employ uj^Jj further to connect 
the pair of Hamiltonians and H2 = H^—vi r , in which 
H?, with edges inside the block, is simply a Hamilto¬ 
nian achieved via translating Hl by L/2 sites. An addi¬ 
tional condition to simulate the effect of vi r here is that 
the effective range of U^Jj should be less than L, i.e., 

^adj should commute with v r i ■ The second form of U a dj, 
U™ = is then yielded and the contained 

|are just eigenstates of Hy . By the same token, 
as the compound systems H^l and H% L are considered, 
the effects of interblock couplings V 12 and V 21 on their 
low-energy behavior could be simulated by imposing the 
transformations, either or £7^., on corresponding 
adjacent segments of the disconnected system H^ + H\ 2 }. 

The analysis above reveals that the regularized basis 
sets of equations III and j2j), or more accustomedly, their 
symmetric and anti-symmetric combinations I'I'ji) = 
\^n') ± I Vfj) and |^ ± ) = |TC,) ± specify ex¬ 


actly the low-energy solutions for compound blocks H 2 L 
and H%l as L —> 00 . It thus illuminates a regularized ver¬ 
sion of the NRG scheme by incorporating the described 
transformation into the algorithm in order to reconstruct 
Hamiltonians for doubly increasing blocks. Historically, 
the speculation through applying a variety of bound¬ 
ary conditions to perform the NRG procedure has ever 
been endeavored by White and Noack in dealing with the 
tight-binding model and the localization problem of one 
dimension . 7], [14]. The discovery of the distinct prescrip¬ 
tion here should enable us to exploit the NRG method 
for general quantum lattice systems, i.e., for interacting 
systems and for systems with high spatial dimensions. 

The key of the present algorithm rests on recognizing 
correspondences among basis states of the Hamiltonians 


with different boundary configurations so as to build the 
transformation U a dj ■ For this purpose, one should sepa¬ 
rate the basis states by quantum numbers and then dis¬ 
cern the matching of states among different sets but with 
identical quantum numbers through comparing their fi¬ 
delities. Due to finite-size effects, the matching of states 
would become less evident as the energy level increases; it 
requests consequently that the initial block to be exactly 
diagonalized should be of considerable size. It should be 
noted that, in spite of the correspondences having been 
established, U a dj’s are still not uniquely determined ow¬ 
ing to the phase uncertainty of the eigenstates. This 
uncertainty, however, doesn’t affect the basic efficiency 
of the scheme since different choices of the phases have 
no influence on the transformations under L —> 00 . In 
practical calculations it can be simply removed, e.g., by 
setting (V’fcl^fc) (or {ip®\ V’D) to be rea ^ anc ^ positive [ 151 ] . 

To verify the efficiency of this new version of the NRG 
scheme, the algorithm has been implemented to calculate 
a periodic spin-1 antiferromagnetic Heisenberg model 
H = Si-S i+ i with 16 sites. The initial Hl and hi + h r 
with L = 8 are exactly diagonalized and the transforma¬ 
tion U^Jj is employed to generate the basis states for 
H%l- The algorithm recovers the exact energy ESE3 to 
at least 5 digits (with to = 47 or more states to be kept). 
As the derived low-lying spectra are shown in Table I, 
some key points to perform the algorithm are recited as 
below. (1) It is convenient to invoke four-component ten¬ 
sors to represent the states: |Vb)® IVy) — where 

the indices ji , j r , j[ and j' r account for the four half blocks 
with size L = 4, respectively. Transformation of U ar ] 0 on 
adjacent halves is then realized by summations over in¬ 
dices ( j r ,j'i ) or (j' r ,ji). (2) Symmetries related to the 
total spin S 2 and S z can be exploited to reduce the com¬ 
puting cost since these invariants are preserved under 
the transformation U a dj- States with definite quantum 
numbers (s, s z ) for the L = 16 block are readily con- 


TABLE I: Energies of 21 low-lying states of the L = 16 
Heisenberg chain obtained by the regularized NRG scheme 
using Presented also includes the information of the 

degeneracy /, quantum numbers of the total spin s and the 
momentum k of the corresponding energy levels yielded by 
the symmetry-preserved algorithm. An accuracy about 1CP 5 
is obtained by keeping m = 47 states (cf. the values of the 
ground and the first-excited energies achieved by exact diag- 
onalization Eo « —22.44681 and E\ ss —22.00401 0,0)- 




s; k 

TO = 21 

to = 47 

Eo (/ 

= 1) 

s — 0; k = 0 

-22.44170 

-22.44639 

Ei (/ 

= 3) 

S = 1; k = 8 

-21.99120 

-22.00326 

E 2 (/ 

= 6) 

s = 1; k = 7,9 

-21.36381 

-21.37456 

e 3 (/ 

= 5) 

s = 2; k = 0 

-21.27190 

-21.28826 

Ei (/ 

= 6) 

s=l;fc= 1,15 

-21.06151 

-21.07610 
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structed by those of L = 8 in virtue of Clebsch-Gordan 
coefficients [l8|, which are then acted on by the trans¬ 
formation Uf d j and are Gram-Schmidt orthonormalized 
afterwards so as to get the desired basis states, g ). 
(3) The algorithm described till now does not involve the 
translational symmetry except for that of the blocked 
translation with 8 sites. To cope with the full transla¬ 
tional symmetry, one should construct states with defi¬ 
nite momentum quantum number k using the basis states 

I ^j,s, Sz ) and = and then di “ 

agonalize the Hamiltonian H% L in symmetry-preserved 
spaces with quantum numbers (s, s z , k). 

Extensions of the proposed NRG scheme to lattices 
with more than one spatial dimension could be naturally 
procured. Take the square lattice as an example: by 
sketching out four different boundary configurations of 
an L x L lattice [see Fig. 1 (a) and (b)], the transforma¬ 
tions U^Jj and U^j could be constructed from eigenstates 
of pair configurations, either (i) and (ii) of Fig. 1 (a) or 
(iii) and (iv) of Fig. 1 (b), respectively. Since the lattice 
Hamiltonian usually possesses symmetries related to a 
specific crystallographic point group, classifying degener¬ 
ate eigenstates via irreducible group representations [19] 
should be a task prerequisite to the procedure. The com¬ 
pound system of the quadrupled block needs to be repre¬ 
sented by a 16-component tensor with indices illustrated 
in Fig. 1 (c). Basis states for the block Hamiltonian with 
periodic boundary conditions are obtained by imposing 
Uadj on each joining region with fourfold adjacent edges, 
which can be realized via summations over corresponding 
indices, respectively. 

As a test case, the two-dimensional scheme is examined 
by a tight-binding model, a single particle hopping on a 
square lattice with a Hamiltonian 

H = ^l(u a l+i,y a x,y + h - c -) + (.vaiy +1 a x ,y + h .c.)l (3) 
%,y 

where a' x y and a XtV are the creation and annihilation op¬ 
erators on the site (x, y) and the hopping parameters are 
set as u = — 1 and v = —3/2. The Hamiltonians of an 
L = 64 square lattice with various boundary configura¬ 
tions are exactly diagonalized of which the eigenstates of 
configurations (ii) and (iv) are recorded as | ipj) and | ip?), 
respectively. Unitary U a dj 's can be built for the present 
system through full correspondences of the derived bases, 
which enables us to calculate states for the quadrupled 
lattice at any energy level. Starting from m states {|^-)} 
closest to a particular energy, the 4m basis states of the 
compound periodic lattice are generated and divided into 
four sets, {|4'~ b ' ± ), j = 1, • • •, in}, in which “ ± ” denotes 
the symmetric or anti-symmetric property of states un¬ 
der t x (y-j = translations along the x ( y ) direction 

by L sites. The lowest energies of the L = 128 system 
achieved by the algorithm via are shown in Table II, 
with which the exact energies and the energies achieved 



FIG. 1: Schematic of configurations of the square lattice, (a) 
Configurations (i) and (ii) of the L x L lattice with open 
boundaries, in which inside edges are either disconnected [con¬ 
figuration (i)] or connected [configuration (ii)]. (b) Lattice 

configurations (iii) and (iv) with connected boundaries. The 
one with disconnected inner edges [configuration (iii)] is equiv¬ 
alent to the configuration (ii) up to a translation along both 
x and y directions by L/2 sites, (c) Configuration of the rep¬ 
resentative tensor for the compound quadrupled lattice. The 
joining regions of adjacent segments are indicated by indices 
{4, 7,10,13}, {8, 3,14, 9}, {12,15, 2, 5}, and {16,11,6,1}, on 
which U a dj should be imposed separately. 


by the primitive NRG are presented for reference. 

Performance of the algorithm employing the prescrip- 
( 2 ) 

tion U ad:l on the above two models, the spin-1 Heisen¬ 
berg chain and the tight-binding model, displays that it is 
slightly less accurate than that utilizing . For exam- 

pie, the algorithm using yields a ground-state energy 
Eq k, —22.43222 for the 16-site Heisenberg chain pro¬ 
vided that the same amount of 47 states are kept. Never¬ 
theless, it is of interest to mention an exceptional case of 
the tight-binding model, where the numerical calculation 
discloses that the symmetric basis states gen- 

erated via imposing Uv.’s on the state vectors of direct 
sum | ifij) ® | ipj) ® | ipj) are exact eigenstates of the 

quadrupled system with periodic boundary conditions. 
In view that the intercepted state vector of any joining 
segments with fourfold adjacent edges represents a basis 
state of the lattice configuration (iii), the transformation 
Uadj 0n each of those joining segments maps it precisely 
onto \4>j), a basis state of the lattice configuration (iv). 
That is to say, direct combinations of basis states of the 
periodic L x L block, ® \tf^) ® \ifij) ® | ip*]), already 
give rise to one fourth of the full set of exact solutions 
of the compound 2 L x 2 L system. This result is indeed 
general: it is independent of the size L and applies also 
to the one-dimensional tight-binding model wherein it is 
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TABLE II: Lowest energies of the L = 128 square lattice with periodic boundary conditions achieved by the regularized 
NRG algorithm via the prescription U^. The energies are calculated in separated spaces of the four sets of basis states 
)ij = 1, • • •, m} with m = 64 and are listed according to their values. The results recover the exact energies to 4-5 digits 
and exhibit correct order, while the results obtained by the primitive NRG scheme are quite inaccurate and fail to yield the 
same ordering. 



{t X ,ty} 

Exact 

Regularized RG 

Primitive RG 


\tx-, ty} 

Exact 

Regularized RG 

Primitive RG 

Eo 

{ + ! +} 

-5.00000 

-4.99960 

-4.99462 

Eg 


-4.99398 

-4.99386 

-4.97853 

E 1 

{“ J +} 

-4.99759 

-4.99732 

-4.99431 

£9 

{+. +} 

-4.99037 

-4.99011 

-4.98698 

e 2 

{~1 +} 

-4.99759 

-4.99728 

-4.98821 

£10 

{+. +} 

-4.99037 

-4.99009 

-4.98314 

e 3 

{ + > _ } 

-4.99639 

-4.99623 

-4.99417 

£11 

{+.-} 

-4.98676 

-4.98674 

-4.98494 

£ 4 

{ + > _ } 

-4.99639 

-4.99617 

-4.98653 

£12 

{+. -} 

-4.98676 

-4.98672 

-4.97730 

£5 


-4.99398 

-4.99395 

-4.99387 

£13 

{+.-} 

-4.98676 

-4.98668 

-4.97698 

£e 


-4.99398 

-4.99391 

-4.98776 

£14 

{ + > — } 

-4.98676 

-4.98666 

-4.96774 

£7 


-4.99398 

-4.99389 

-4.98463 

£15 

{+. +} 

-4.98555 

-4.98541 

-4.97742 


very easy to verify. 


In the present algorithm at least two truncated basis 
sets are kept and they should be expressed in a com¬ 
plete set of spin bases. This should be the case at each 
iteration since any such truncated set is too incomplete 
to represent other basis states with different boundary 
configurations. As a result, repetition of the described 
RG procedure is constrained by the capacity of storing 
vectors with exponentially growing dimensions. In this 
sense, the regularized NRG scheme is applicable to finite 
lattices rather than the infinite one since the iterative 
performance of the procedure for the latter system will 
break down inevitably. This indeed offers the underly¬ 
ing cause which supports partly the folklore stating that 
all real-space RG schemes are necessarily inaccurate for 
lattice systems. 


Alternatively, if one chooses to sacrifice some of the ac¬ 
curacy, an iterative procedure could be achieved, at least 
for one-dimensional systems, by projecting the derived 
basis states, e.g., those of H 2 l and H% L , on the tensor 
product space of the to states kept to get bases \ipj) and 
|i// J ) of dimension mr. The prescribed transformation, 


U a dj , could be built by virtue of these basis states which 
is then imposed on the states | ipj) ® | ipj') to generate ba¬ 
sis states for the pair of Hamiltonians with larger size. 
The truncation of the algorithm here is distinctly differ¬ 
ent from that of the primitive NRG scheme in view that 
the to states kept are taken from a Hilbert space of di¬ 
mension to 4 instead of dimension m 2 . It is also worthy to 
note that the strategy adopted here, projecting the ba¬ 
sis states of the compound block onto a restricted space 
of tensor product states of subsystems, is what has been 
done in the contractor RG method 20H22]. In compari¬ 
son, the interblock coupling of the effective Hamiltonian 
in the latter scheme is obtained by substracting the con¬ 
tributions of contained subclusters, while in the present 
scheme the effective Hamiltonian of the compound sys¬ 


tem could be constructed directly by virtue of the regu¬ 
larized basis states. 
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